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Abstract. 

Detection and analysis of Wolf-Rayet (WR) star populations in young 
starbursts is a powerful tool to reveal the properties of the star-forming 
t^J- , events. We present new high S/N optical spectra of a sample of well- 

ed ' known WR galaxies, in which we detected for the first time the spectral 

signature of WC stars. The detailed comparison of WN and WC star 
populations derived in these galaxies, with the most recent model pre- 
dictions, allows to put tight constraints on starburst parameters, such as 
' the age and duration of the burts, and the IMF. 
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Introduction 



Starbursts play a major role in the global process of galaxy formation and evo- 
^ . lution. A better knowledge of their intrinsic properties (duration, IMF, etc) in 

5-j | nearby starburst galaxies is a key issue to follow the evolution of the large pop- 

ulation of star-forming galaxies observed up to redshift z ~ 4 (e.g. Steidel et al. 
1998). The so-called "Wolf-Rayet (WR) galaxies" may be the ideal laboratories 
to derive starburst properties since these objects harbour the most massive stars 
known, O stars and their descendents WR stars, which allow to probe the upper 
part of the IMF and the youngest stellar populations. 

Since the initial compilation of Conti (1991) listing 37 galaxies, the number 
of known WR galaxies has grown rapidly to about 130 in the last compilation 
of Schaerer, Contini &; Pindao (1999, hereafter SCP99). While most of them 
fall in the category of "classical" star-forming galaxies (Hn galaxies or Starburst 
Nucleus Galaxies), broad WR emission lines have also been detected in more 
"exotic" objects like Seyfert 2 galaxies (e.g. Heckman et al. 1997; Kunth & 
Contini, these proceedings) or giant cD galaxies, located in the centre of galaxy 
clusters with strong cooling flows (Allen 1995). 

The presence of a large number of WR stars has already provided interesting 
informations about star formation in WR galaxies, indicating recent (<^ 10 Myr) 
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star-forming events and the existence of massive stars (e.g. Vacca & Conti 1992, 
hereafter VC92). But in practice, the detection of the broad emission line Hen 
A4686, attributed mainly to WR stars of WN subtype, was often simply used. 
The strongest emission line of WC stars, C iv A5808, has quite rarely been 
observed so far in WR galaxies. Although Conti (1991) claimed the absence of 
convincing evidence for carbon features in WR galaxies, there is no doubt about 
their existence anymore since about 56 objects showing both broad Hen A4686 
and Civ A5808 features are listed in SCP99. 

The existence of WC stars in WR galaxies, although possibly less numerous 
than WN stars, is indeed expected both from observations of WR populations in 
the Local Group (e.g. Massey & Johnson 1998, hereafter MJ98), and from stellar 
evolution models (Maeder & Meynet 1994). Surprisingly, the homogeneous, 
fairly high S/N spectra from the sample of VC92 containing 12 regions with 
Hen A4686 detections revealed only one Civ A5808 detection (in He 2-10 A). If 
true, this would certainly contradict the predictions from synthesis models for 
young starbursts (Meynet 1995, Schaerer & Vacca 1998, hereafter SV98), which 
otherwise compare well to observations in the Local Group (Maeder & Meynet 
1994). To verify if this apparent discrepancy really holds we have searched 
for WC stars in well-known WR galaxies (NGC 3049, He 2-10, NGC 3125, 
NGC 5253, Tol 89). Previous accounts of our results are given in Schaerer et al. 
(1997) and Schaerer, Contini & Kunth (1999, hereafter SCK99). 

2. Observations and spatial analysis 

Long-slit spectra of galaxies were obtained on the nights of 1995 April 24 - 26 at 
the ESO 2.2m telescope. The data were acquired with the EFOSC2 spectrograph 
which gives a spectral coverage of 4400 - 6500 A with a resolution of ~ 5 A. 
The slit was oriented in order to cover regions with previous detections of WR 
stars. The total integration time for each galaxy ranges from 80 to 120 min. 
The seeing was relatively stable during the observations with a mean spatial 
resolution of about 1". The procedure for data reduction is fully described in 
SCK99. 

We used our long-slit spectroscopic observations to compare, for each galaxy, 
the spatial distributions of nebular and WR emission-line intensities and the 
stellar continuum emission. The results for each galaxy are given in SCK99; 
they can be summarized as follows. There is an offset between the distribution 
of nebular lines and stellar continuum in He 2-10 and Tol 89. The WR emis- 
sion features follow the stellar continuum with no significant shift between He n 
A4686 and C iv A5808. The spatial extension of WR features ranges from ~ 
60 pc (NGC 5253) to ~ 270 pc (He 2-10 and NGC 5253). Two distincts WR 
regions are identified in NGC 3125 and NGC 5253. 

3. The WR regions 

One-dimensional spectra were extracted corresponding to the different regions 
where spectral signatures from WR stars have been detected. Normalized spec- 
tra showing the region around the "blue" and "red" WR bumps are given in 
Fig. [l]. Several broad emission lines due to WR stars have been detected over 
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Figure 1. Spectra of galaxies showing the "blue" (left) and "red" 
(right) WR bumps 

the wavelength range covered by our spectra: the blend of N in A4640 and C in 
A4650, Hen A4686, and Civ A5808. Equivalent widths (EWs) of these emission 
lines are listed in Table [l], together with the basic nebular parameters (metal- 
licity, extinction, etc). 

One of the main results of these observations is the unambiguous detection 
of broad (FWHM ~ 50 - 90 A) C iv A5808 emission in all WR regions of the 
observed galaxies (see right panel of Fig. [I]), which clearly indicates the presence 
of WC stars in these regions. Whereas all the objects were previously known 
to show the broad WR bump around 4700 A, our spectra are the first to show 
the red WR bump around 5800 A. As already noted (Schaerer et al. 1997, 
SCK99), a sufficiently high S/N (^ 40) is required for the detection of such a 
C iv A5808 line in integrated spectra of WR galaxies. This explains why this 
line has remained undetected in previous ground-based spectroscopical studies, 
except for He 2-10 (VC92). 

In all our seven emission-line regions we thus find signatures of both WN 
and WC stars. The dominant WR subtypes in each region can be constrained 
using measured line widths and intensity ratios between WR emission lines (see 
SCK99 for a full description). All regions (except the high-metallicity galaxy 
NGC 3049) contain a mixed population of WNL, and early WC and/or W03-4 
stars. In NGC 3049, the detection of Cm A5696 indicates the predominance 
of WCL stars. The presence of intermediate-type WN/WC stars cannot be 
excluded. 
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Table 1. The broad WR emission lines 



Line / Galaxy 


He 2-10 NGC 3049 


NGC 3125 


NGC 5253 


Tol 89 




A 


B 


A 


B 




EW(np) [A] 


23± 3 44± 3 


93± 9 


70± 7 


269± 23 


112± 10 


68± 17 


(O/H) [Z ] 


0.42 1.20 


0.17 


0.20 


0.20 


0.20 


0.25 


CH/3 


0.56 0.23 


0.40 


0.64 


0.44 


0.20 


0.18 


EW [AJ OF BROAD EMISSION LINES FROM WOLF-RAYET STARS 


N in A4640 


2.7±0.5 5.2±0.5 


3.2±0.5 


2.8±1.0 


2.1±1.3 


1.0±0.5 


3.4±0.5 


+Cm A4650 














Hen A4686 


2.3±0.5 5.6±0.5 


7.3±1.0 


2.8±0.5 


3.0±0.7 


3.3±1.1 


6.1±0.5 


Cm A5696 


2.0±0.5 












C iv A5808 


1.5±0.5 4.7±0.6 


3.1±0.4 


6.2±1.0 


2.5±0.6 


6.0±1.1 


11.9±1.5 



4. Populations of WR and O stars 

The approximate number of WN, WC, and O stars derived from the observed 
emission lines are listed in Table |2[ The number of WN and WC stars is calcu- 
lated from the luminosity of the WR emission lines, assuming that the dominant 
contributors to the broad He 11 A4686 and C iv A5808 lines are respectively WNL 
and WC4 stars (see previous section). Under the condition of case B recombi- 
nation and assuming that all the ionizing photons emitted by the stars are 
absorbed by the gas, the total number of Lyman photons can be derived from 
the observed luminosity of the H/3 emission line. The number of O stars is thus 
estimated, taking into account the ionizing photon contribution from WR stars, 
the age of the stellar population and the IMF (see SCK99). The number of O 
stars found for the WR regions are between ~ 500 and 7000. The case of NCG 
3049 is somewhat exceptional (see SCK99 for a discussion). According to the 
models, the bulk of the ionisation in this galaxy is provided by late O and/or 
B type stars and WR stars. Uncertainties on absolute number of WR and O 
stars reported in Table ^ are quite large. The main sources are the following 
(see SCK99 for details). There is a large uncertainty on the average luminosity 
of WNL stars in the He 11 A4686 line, and to a lesser extent of WCE stars in the 
C iv A5808 line. The dominant WC subtype is not well constrained. There may 
be a nebular contribution to the He 11 A4686 emission in NGC 3125-A and NGC 
5253-A (see Fig |l]). A significant fraction of ionizing photons may escape from 
Hll regions. 

The number ratios WR/O and WC/WN are reported in Table pi We obtain 
WR/O number ratios of ~ 0.03 — 0.4 (and even larger for NGC 3049), systemat- 
ically higher than the predictions for constant star formation at the appropriate 
metallicity (Maeder & Meynet 1994), but within the range of instantaneous 
burst models with different IMF slopes (SV98). A trend of increasing WR/O 
ratios towards higher metallicity is found as expected (e.g. Meynet 1995). The 
number ratio of WC/WN stars is typically 0.2 — 0.4, with no systematic varia- 
tion with metallicity. For the low- metallicity objects {Z ~ 0.2 Z Q ; NGC 3125, 
NGC 5253, and Tol 89) the derived WC/WN ratio is larger than what is found 
in the Local Group (except IC 10) at similar metallicity (MJ98). At higher 
metallicities, He 2-10 and NGC 3049 show, however, WC/WN ratios below the 
trend observed by MJ98. 
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The finding of a fairly constant WC /WN ratio may seem surprising at first 
sight. Contrary to regions of constant star formation, likely representative of the 
Local Group samples, regions of short star formation could a priori show quite 
a large range in WC/WN, depending on the age of the starburst. However, as 
explained in SCK99, burst durations of ^ 2 Myr already suffice to smooth out 
the WC/WN ratio. This may well explain the small range of WC/WN found 
for the low-metallicity objects. More surprising is the low WC/WN ratio of the 
remaining higher metallicity objects (He 2-10 and NGC 3049) compared to the 
WC/WN value of MJ98 at a similar metallicity. Indeed, the probability of find- 
ing WC/WN ratios below the equilibrium value attained in regions of constant 
star formation should be quite small. However, if the number of WC stars is 
systematically underestimated by a factor of 3, our observations may well all be 
larger than the observed WC/WN trend with Z. We conclude that the WC/WN 
ratios of the low-metallicity galaxies can be understood quantitatively with burst 
models of reasonably short but non-zero duration. Additional observations of 
WC and WN populations, especially for regions of higher metallicities, would 
be very helpful. 



Table 2. Populations of WR and O stars in WR regions 



Galaxy 




Age burst 
[Myr] 


WN 


WC 


O 


WR/O 


WC/WN 


He 2-10 




5.5 - 6.0 


1100±520 


> 250 


2450 - 4900 


0.37 


0.23 


NGC 3049 




5.5 


510±240 


> 170 






0.33 


NGC 3125 


A 


4.5 - 5.0 


500±230 


> 70 


3240 - 6470 


0.12 


0.14 




B 


4.5 - 5.0 


530±250 


> 200 


3450 - 6900 


0.14 


0.38 


NGC 5253 


A 


3.0 


26±13 


> 9 


960 - 1080 


0.03 


0.35 




B 


5.0 


27±13 


> 10 


680 


0.05 


0.37 


Tol 89 




4.5 - 5.0 


240±110 


> 150 


640 - 1270 


0.41 


0.62 



5. Age, duration and IMF of starbursts 

One can try to constrain the main burst parameters (age, IMF, star formation 
history) by comparing the observed WR signatures to the predictions from the 
latest population synthesis models for young starbursts (SV98). We consider 
two basic free model parameters: the IMF slope and the duration of the star- 
forming event. The metallicity Z and the starburst age (from EW(rlf3)) are 
adopted from the observations. A power law with a slope a is adopted for the 
IMF (in our notation a = 2.35 for a Salpeter IMF). The upper mass cut-off 
is generally set to M up = 120 M . We consider burst models with different 
durations At starting at time t = 0, with the limiting case of an instantaneous 
burst (At = 0). SV98 models are based on the latest Geneva stellar evolution 
tracks and only single-star models are considered. A detailed discussion on the 
assumptions or uncertainties of the models is given in SCK99. 

The comparison with the observations are performed in two ways. One 
both compare relative line intensities of the WR features with respect to the 
nebular H/3 emission (in short WR/H/3), and equivalent widths of the WR lines. 
Observationaly the determination of these quantities can mainly be "perturbed" 
by three effects: 
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(i) We do not count all the ionizing photons produced by massive stars. This 
may happen if the nebular emission from the Hll region of interest is not 
entirely included in the slit, and/or if ionizing photons can escape from 
the Hll region. 

(ii) Stars and gas suffer from a different extinction. 

(Hi) An underlying older population contributes additional continuum light. 

The relative intensities WR/H/3 are affected only by (%) and (ii). Effect (i) in- 
creases the relative WR/H/3 intensities. The same holds for (ii) if the stellar 
light is less extincted than the gas (cf. e.g. Calzetti 1997, Mas-Hesse & Kunth 
1999). The EWs are only affected by (Hi), which decreases the observed value. 
We refer to SCK99 for a discussion on the quantitative importance of these 
effects. Given these potential difficulties both comparisons of relative line inten- 
sities WR/H/3 and WR equivalent widths have been performed for all objects. 
Significant differences between the two methods likely indicate some difficulty 
with (i) and (Hi). In this case the simplest meaningful comparison is between 
the observed WR equivalent widths and the maximum value predicted by the 
models (irrespectively of the age); it only depends on (Hi). 

In Fig. |^ we compare, for all the WR regions, the observed relative inten- 
sities and EWs of WR lines with the model predictions of SV98 at metallicities 
Z = 0.02 (solar), 0.008, and 0.004 for an instantaneous burst with a Salpeter 
IMF. In both plots, W(H./3) is used as a time indicator (decreasing W with 
time) in the synthesis models. Comparison of the individual objects (classified 
by metallicity) with the model predictions can be summarized as follows. 

5.1. High-metallicity galaxies 

The predictions of the standard model (instantaneous burst, Salpeter IMF) are 
in agreement with the observed WR line strengths of the high-metallicity galax- 
ies, particularly for NGC 3049. For He 2-10, the observed EWs are considerably 
smaller than the maximum value predicted by the models, but the relative line 
intensities are close to the predicted maximum. This difference could be at- 
tributed to the shift between gaseous and stellar emission (see sect. ^]). 

5.2. Low-metallicity galaxies 

For the low-metallicity objects, the comparison between data and model pre- 
dictions is much more complex. Only some of the observed WR line intensities 
and EWs lie well on the predicted model curves plotted in Fig. ||. For the con- 
sidered objects we dispose in total of 15 measurements of broad emission lines 
(all lines and all regions). The relative line intensity of 7 WR regions is above 
the maximum value predicted by the standard model. On the other hand the 
majority (13 of 15) of the EWs are well in the range of the model predictions. 
The two remaining cases are VK(Heii A4686) of NGC 3125-A and Tol 89, which 
could well be overestimated due to nebular contamination and contributions 
from other lines in the blended region, and an exceptionally large C IV A5808 in 
Tol 89. The uncertainties affecting the observed quantities are fully discussed in 
SCK99. Only for Tol 89 we have direct evidence that the displacement between 
gas and stars is likely of importance (see sect. ^]). For this galaxy we measure 
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Figure 2. Observed vs. predicted WR line intensities relative to H/3 
(left) and WR equivalent widths (right) 



an extraordinarily strong C iv A5808 emission (see Table [T]) which exceeds the 
values shown in Fig. [2|. We have no simple explanation for this strong emission. 
Other observations are required to confirm the measurements for Tol 89. 

One possibility to explain the observations of some regions with large WR/H/3 
intensity ratios and large W^(WR) may e.g. be a flatter IMF. Models with differ- 
ent burst durations and IMF have been tested in SCK99 (their fig. 12). These 
show that, within reasonable assumptions for the IMF slope, the duration of 
star formation episode is limited to At <^ 2 — 4 Myr. What can be said about 
the IMF slope? As mentioned earlier the observations of the WR features in 
NGC 3049, He 2-10, and NGC 5253 are quite compatible with a Salpeter IMF. 
However, the large observed values of 1^(4686) in NGC 3125-A and Tol 89, and 
VF(5808) in Tol 89 may require more WR stars. This can be obtained by in- 
voking a flatter IMF, which increases the WR equivalent widths. However, we 
note that intrinsicly the model predictions for He n A4686 are not very sensitive 
to changes of the IMF slope. More IMF-sensitive are the features of WC stars, 
which descend on the average from more massive predecessors than WN stars 
and from a narrower range of initial masses. Although the standard models show 
some deficiencies, it is difficult to claim significant differences with respect to a 
Salpeter IMF (see SCK99 for details). We therefore conclude that within the 
uncertainties our observations are compatible with a Salpeter IMF. Although 
we cannot exclude this possibility, no clear case requiring a significantly flatter 
IMF is found. Much steeper IMF slopes may, however, not be compatible with 
our data. 



6. Conclusions 

The main result of this work is the unambiguous detection of WC stars in five 
previously known WR galaxies, which brings to ~ 54 the total number of ex- 
tragalactic objects known to harbour both WN and WC stars (cf. SCP99). The 
relative weakness of Civ A5808 compared to Hen A4686, and its larger width 
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requires sufficiently high S/N (;> 40) to be detected. This explains most likely 
the non-detection in previous observations. 

From the observed WR emission lines we concluded that all WR regions 
(except NGC 3049) contain a mixed population of WNL, and early WC and/or 
W03-4 stars. Weak Cm A5696 indicative of late WC stars is found in the 
high-met allicity WR region of NGC 3049. The estimated WC/WN number 
ratios (typically between 0.2 - 0.4) show no clear trend with metallicity. For 
low-met allicity objects (Z ~ 1/5 Z Q ), these values are larger than the observed 
WC/WN ratios in Local Group objects with similar Z (MJ98). We argue that 
our WC/WN values are compatible with expectations for regions of short star 
formation. For He 2-10 and NGC 3049 the derived WC/WN ratio is below the 
trend given by MJ98. This can be due to the fact that we observe a short 
burst at a particular time (quite low probability), or that the number of WC 
stars is systematically underestimated. The solution awaits new observations 
and quantitative analysis of "WC+WN galaxies" of different metallicities. 

Detecting both WN and WC features provides a considerable improvement 
to derive burst parameters. Given the WN — > WC sequence, which is expected 
to be followed only by the most massive WN stars, it is clear that the predictions 
of WC/WN populations are particularly sensitive to the evolutionary scenario 
and burst parameters (e.g. IMF, burst duration). 

The comparisons of the observed EWs and line intensities relative to H/3 
with models do not all show a simple picture. The most important effects and 
uncertainties which may affect such a comparison have been amply discussed 
in SCK99. In order to reproduce the WR lines, young populations with short 
durations of star formation are required. From our quantitative modeling, we 
find a conservative limit for burst durations of typically At <^ 2 — 4 Myr. The 
majority of the observed quantities can be reproduced reasonably well by the 
SV98 models with a Salpeter IMF. Although some WR lines may indicate a 
flatter IMF in some regions, no clear case requiring a significantly flatter IMF 
is found. Much steeper IMF slopes may, however, not be compatible with our 
data. These results are in agreement with other studies of similar objects (see 
SCK99 and references therein). 
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